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I. INTRODUCTION 



Since Wilson's pioneering work there has been considerable activity to solve field 
theories non-perturbatively by numerical simulation on Euclidean spacetime lattices. Today 
lattice gauge theory is a mature field, and current state-of-the-art lattice QCD calculations 
are beginning to confront the challenges provided by the hadron spectrum. For an overview 
of lattice methods, see [2J. One aspect to these numerical simulations is the finite-size scaling 
of observables. The finite spacetime volume employed on the lattice is a source of systematic 
error in the numerical determination of observables. Thus the study of field theories in finite 
volume, while a theoretical curiosity, is also of practical utility. 

Recent work [3[ suggests that electromagnetically gauge invariant amplitudes at finite 
volume may differ from their infinite volume form. Specifically investigated was the finite-size 
scaling of nucleon electromagnetic and spin polarizabilities that arise in nucleon Compton 
scattering (see, e.g., j3, @]). A goal in ^ was to address systematic errors in the extraction of 
polarizabilities from classical b ackg round field methods employed in lattice simulations 0, 
0, H, M, US, EH El, 13, 14, 15, 16j . An analysis of the finite volume behavior of nucleon 



polarizabilities was presented, as was an oddity relating to the zero-frequency scattering 
amplitude. In infinite volume, the zero-frequency Compton amplitude is fixed by gauge 
invariance to be proportional to the total charge squared. Finite volume modifications, 
however, were found for nucleon Compton scattering at zero frequencies [3]. These results 
suggest a finite volume renormalization of the basic interaction between the photon and the 
hadron's charge. In this work, we show that gauge invariance in finite volume allows for 
such modifications to zero-frequency photon couplings. In essence, conserved currents are 
not protected from additive renormalization as they are in infinite volume. For definiteness, 
we focus on the chiral dynamics of pions coupled to photons [13], but could just as well 
choose any interacting field theory coupled to gauge fields. 1 

Our presentation is organized as follows. First in Sec. [Til we analyze the electromagnetic 
interactions of pions in finite volume. We demonstrate the infrared running of electromag- 
netic current matrix elements by explicit one-loop calculations in chiral perturbation theory 
(xPT). In Sec. 1111} consequences of gauge invariance on a torus are detailed. Gauge invariant 
zero-mode interactions allow for infrared renormalization of electromagnetic couplings. We 
write down gauge invariant, zero-frequency effective field theories for pions that reproduce 
our one-loop finite volume %PT results. Understanding such volume effects is necessary in 
practice for the extraction of infinite volume physics from lattice QCD simulations. We show 
how our results are consistent with Ward identities and low-energy theorems in Sec. IIVI A 
conclusion in Sec. |V] summarizes our findings, while a glossary of finite volume functions is 
provided in the Appendix. 



II. PIONS IN FINITE VOLUME 

The chiral Lagrangian is written in terms of a coset field £ = exp(2i$//) which 
parametrizes the Goldstone manifold arising from spontaneous chiral symmetry breaking: 



An instructive alternate example is the QED electron. Straightforward evaluation shows that the electron 
vertex function at zero frequency is modified by volume effects. This modification, however, is infrared 
divergent and we have chosen to avoid such difficulties by using a theory that is infrared finite. 
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FIG. 1: One-loop graphs required to evaluate the pion electromagnetic current. On the left 
appears the wavefunction correction; while, on the right, diagrams contributing to the pion form 
factor. Vertices shown are generated from the leading-order xPT Lagrangian, Eq. ([2]). 

SU(2) L ® SU{2) R — > SU{2) V - The pions are contained in the matrix $, explicitly as 



In our conventions, the dimensionful parameter / = 132 MeV. The chiral Lagrangian provides 
an effective theory of low-energy QCD. At leading-order in an expansion in momentum, p 2 , 
and quark mass, m q , there are two terms in this Lagrangian 



where m q is the quark mass matrix, m q = diag(m u ,mA. We shall work exclusively in the 
isospin limit, m u = rrid = m. The kinetic term of the chiral Lagrangian includes a U(l) 
gauge covariant derivative that couples pions to photons, D^T, = c^E + ieA^ where 
the quark electric charge matrix, Q, is given by Q = diag(2/3, —1/3). 

Expanding the Lagrangian in Eq. to tree level, one sees that the pions are correctly 
normalized and their mass, m n , is given by = 2Xm. The couplings of pions to zero- 
momentum photons at tree level can be read off from Eq. (J2J), from which we find their 
canonical charges. We now investigate whether loop corrections in a finite spatial volume 
modify these couplings. 

A. Charged pion current 

To consider the one-loop corrections to the electromagnetic current of charged pions, we 
accordingly expand the x?T Lagrangian in Eq. (j2J) to second order to generate vertices 
for one-loop graphs. Furthermore local terms at higher-order can then contribute at tree- 
level, but these are absent for zero-frequency photons. Thus we need to determine only the 
diagrams depicted in Fig. [fl 

In the limit of zero frequency and infinite volume, the current matrix element between 
charged pion states is required by gauge and Lorentz invariance to be 




(1) 




(2) 




(3) 
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where the overall sign reflects the charge of the pion. It is a straightforward exercise to verify 
the above form at one- loop order in infinite volume. At an intermediate step, we reach the 
result 

from which we see the wavefunction correction to the tree-level vertex is exactly canceled 
by the loop contributions to the form factor at zero frequency [HI]. In this way, the matter 
fields do not contribute to the running of the coupling and Eq. ([3]) is preserved. 

In finite volume, we repeat the calculation of the pion current to one-loop order. We 
consider each of the three spatial directions of finite length L, and the quark fields subjected 
to boundary conditions that maintain discrete translational invariance. For definiteness, we 
assume periodic boundary conditions. 2 As the pions are point-like objects in the effective 
theory, they satisfy the same boundary conditions as the point-like interpolating field ~ 
q(x)j5q(x). Pions are hence also periodic with quantized spatial momentum modes of the 
form 

q = —n, (5) 

where n is a triplet of integers. To keep matters simple, we keep the temporal extent 
infinite as is commonly done to determine finite-size effects for lattice QCD observables. 3 
To evaluate the pion current, we use the finite volume theory defined by Eq. (TSJ). The loop 
diagrams shown in Fig. [1] are again generated. The only difference compared to infinite 
volume is that the spatial momenta of real and virtual states are quantized. The finite and 
infinite volume theories share exactly the same ultraviolet divergences, so we can calculate 
the finite volume effect by matching the two theories in the infrared. For an observable X 
calculated in both finite, X(L), and infinite, X(oo), volumes, we have 

X(L) = X(oo) + AJf(L), (6) 

where the matching term, AX(L), is free from ultraviolet divergences and gives the finite 
volume effect. 

Returning to Eq. we can carry out the finite volume matching, Eq. (Q, for the pion 
current. We find 

(vr ± (P)|J / ,|7r ± (P)) = ±2eP, jl - $± [2I 1/2 (m*,L) + m 2 / 3/2 (m 2 , L)] | , (7) 

where Ip(m 2 ,L) is defined in the Appendix. Results are consistent with charge conjugation 
invariance and the current is only modified in the spatially finite directions. Specifically the 
virtual pion cloud in finite volume screens the current of the infinite volume pion. In Fig. [2, 



2 Similar results for anti-periodic boundary conditions, for example, can be derived easily using a modified 
momentum quantization condition for the quark fields. As pions remain periodic, the expressions we 
derive also hold for anti-periodic quarks. 



We implicitly choose m T L > 1 so that pion zero modes do not become strongly coupled 19|, l20(. With this 



assumption, the ordinary xPT power counting in infinite volume can be carried over to finite volume 21 1 . 
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FIG. 2: Finite volume screening of the pion current. The relative difference in pion current A J n + 
is plotted as a function of the box size L, for a few values of the pion mass. 

we plot the finite volume modification to the pion current. Here the relative difference in 
the current matrix element at finite volume versus infinite volume, AJ V +, given by 



AJ n+ = 



(n + (P)\e-J 


vr + (P)) L -{ 


7r+(P)|e- J 


7T + (P))oo 


(7r+(P)|e- J 


vr+(P))oo 



(8) 



is plotted as a function of the length L of the spatial dimension. We have used a unit vector 
e to project onto the spatial part of the current. Accordingly the pion cannot be at rest, 
P ^ 0. Subscripts on matrix elements denote the box size, with infinity corresponding to 
infinite volume. The finite volume effect is exponentially suppressed in asymptotic (m n L 3> 
1) volumes. Consequently taking the infrared cutoff, 1/L, to zero, the additive current 
renormalization vanishes and infinite volume limit is maintained. 

B. Neutral pion current matrix elements 

Charge conjugation invariance demands the identical vanishing of single current matrix 
elements between neutral pion states. Indeed whether the calculation of the neutral pion 
current is carried out in infinite or finite volume, we find zero for the matrix element. The 
577(2) flavor structure of the form factor diagrams shown in Fig. [T] ensures this vanishing 
and consistency with charge conjugation. 

Neutral pion matrix elements of an even number of electromagnetic currents, however, 
are not restricted to vanish by charge conjugation invariance. Indeed, it is well known that 
the neutral pion has electric and magnetic polarizabilities that can be predicted at one-loop 
order in %PT solely in terms of / and m n [22]. Such polarizabilities arise at second order in 
the low-frequency expansion of the matrix element of two currents (the so-called Compton 
scattering tensor). The Compton tensor also has a term at zeroth order in the photon 
frequencies 



which is sensitive only to the longest ranged electromagnetic interaction. This term in 
the Compton tensor, when combined with relevant phase space factors, yields the classical 
Thomson scattering cross section, a = 87r(<5e) 4 /3m^. For the neutral pion, the total charge 
is zero and the longest ranged interaction vanishes. 




(9) 
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FIG. 3: One-loop contributions to neutral pion Compton scattering in xPT. Vertices shown are 
generated from the leading-order Lagrangian. 
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FIG. 4: Zero frequency Compton amplitude for the neutral pion. The finite volume amplitude AT 
is plotted as a function of the box size L, for a few values of the pion mass. In infinite volume, this 
amplitude is identically zero. 



Using the Lagrangian denned in Eq. ([2]), we can determine the Compton amplitude 
for pions. We restrict our attention to the zero-frequency amplitude. Due to charge neu- 
trality, there are no tree-level couplings to the neutral pion. At one-loop order, evaluation 
of the diagrams shown in Fig. [3] is required to determine the Compton amplitude. At an 
intermediate step in the calculation, contributions from all six diagrams can be simplified to 

T» v (u -u-0)- —p- J ^ [q2 _ -p; , (10) 

which vanishes. Hence in infinite volume, a delicate cancellation between all diagrams 



maintains the vanishing of the Compton amplitude at zero frequency [23|, [24J. On the other 
hand, the same is not true in finite volume. Carrying out the one-loop matching between 
finite and infinite volume theories, Eq. ([6]), for the Compton amplitude in Eq. ([TO]) , we find 



m 



■i 



T^(cu = J = 0) = e 2 (V^o - g^)^ J 5/2(<, L). (11) 

Thus when one considers the purely spatial components of the Compton tensor, the neutral 
pion has an effective charge- squared, cf. Eq. ([9]). In transverse gauge, the above expression 
gives the amplitude to scatter zero-frequency photons off the neutral pion. There is a non- 
vanishing contribution to this scattering amplitude when the pion is confined to a periodic 
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box with size on the order of the pion Compton wavelength. When the box size becomes large 
compared to this scale, the amplitude is exponentially suppressed and infinite volume results 
are recovered. We demonstrate this in Fig. HI where we plot the finite volume amplitude 
AT defined by 

AT = ~g^T liV (u = u' = 0). (12) 

III. GAUGE INVARIANCE ON A TORUS 

To explain our above results, we investigate electromagnetism in finite volume. The 
analogous finite temperature case is well known and described, e.g., in (25[. Because our 
applications are with classical background fields, or equivalently current operator insertion 
methods in lattice field theory, there are no quantum corrections to the photon field itself. 4 



A. Spatial Torus 

Let us consider a classical electromagnetic field defined on a finite spatial torus with 
infinite time extent. On the gauge field A^(x,t), we impose periodic boundary conditions 
and expand A^x^t) in Fourier modes 

A M (x,t) = ^i M (n,t)e 2 —/ i , (13) 

n 

where n = (n x ,n y ,n z ) e Z 3 . It is convenient to separate out the zero-mode contribution, so 
we write 

A li (x,t) = A li (t)+A ll (x,t), (14) 

where the zero mode A^(t) = A^(0,t). 

Under a gauge transformation, the photon field transforms in the familiar way, 
A^x, t) — > Afj,(x, t) + dfj,a(x, t), and observables are invariant. Requiring the gauge trans- 
formed field to be single valued mandates that d^a^x, t) is periodic. Thus we can decompose 
the gauge function a(x,t) into the sum of two terms, a(x,t) = a (x,t) +a(x,t), where 

a (x, t) = a (t) + ex • x, (15) 

and 

a(x, t) = J2 "( m ' *) e 2mm - x/L . (16) 

Here we have dropped all overall irrelevant constants, and the vector ex is a constant vector. 
Using this decomposition for the gauge function, the gauge field transforms as 

A^t) — ► Ap(t) + d^x, t) ^ 
A^(x,t) — ► A^(x, t) + d^a(x, t) ' 



Dynamical photons in QED cause additional complications as the vector current is renormalized in infinite 
volume 20| . With classical background fields, penguin graphs are absent and such renormalization does 



not occur. 
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In particular, the photon zero-mode transforms as 

M t) _ AM + I*" ' ' " = ?. (18) 
[at, ii — i 

The time-component of the zero mode is absent from the field strength tensor. The remaining 
three components of the zero mode field are translated by a constant under the gauge 
transformation. In the gauge invariant free theory, each spatial component of the zero mode 
is thus a massless one-dimensional scalar. 



B. Coupling to Matter 

For a generic matter field ip(x,t) of unit charge, the effects of a gauge transformation 
show up as a local phase factor 

(p(x,t) — ><p'(x,t) = e~ l ' a(a! '*V0M)- (19) 

Now we assume that the matter field (p(x,t) is subject to periodic boundary conditions. 
We again split the gauge function into zero mode and non-zero mode pieces, a(x,t) = 
ao(x,t) + a(x,t). With the form given in Eq. (fT6|) . we see that the non-zero modes will 
maintain the periodicity of the matter field under the transformation in Eq. ([19]) . The same 
is not in general true of the zero modes given the form of ao(x, t) in Eq. (Tl5l) . If the gauge 
transformed matter field is to remain periodic under translations by L, then we must have 
the quantization condition 

a = — n, (20) 

on the spatial zero mode part of the gauge function. This quantization condition reduces 
the continuous translational invariance of the spatial zero modes to discrete translations. 
As gauge transformations are now less general compared to infinite volume, more gauge 
invariant operators can be built. 

Imagine that we start with some microscopic theory with electromagnetic interactions. 
Take the scalar field <^(x, t) as some composite low-energy degree of freedom of this the- 
ory. Further we assume that the energies of interest are ultra-low in the sense that any 
interactions of <p(x,t) with itself or other fields have been integrated out. In the absence of 
electromagnetism, e.g., we have a simple single particle effective theory 5 

C= \d^\ 2 -m{Lf\^\\ (21) 

where m(L) is a running mass that depends on the infrared cutoff 1/L (and parametrically 
depends on the other couplings, masses, etc. that have been integrated out of the theory). 
Running the cutoff to zero completes the infrared sector of the theory and produces infinite 
volume physics. 



5 We have written only 5*0(4) symmetric terms in Eq. (|2Tj) . Strictly speaking volume corrections will reduce 
the dispersion relation down to only cubic symmetry. For pions in finite volume the first 5*0(4) breaking 
effects occur at two-loop order in the chiral expansion. 
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Now we include electromagnetism in this single particle effective theory by adding all 
possible gauge invariant operators. The minimal coupling prescription, — ■> = d^ + iA^, 
renders the kinetic term of Eq. (T2T]) gauge invariant. Because we imagine ip is a composite 
particle, there can be non-minimal couplings that respect gauge invariance, e.g., the E 2 and 
B 2 terms for the particle's polarizabilities. Further terms in this ultra low-energy theory 
are allowed, however, because SO (4) is not respected, and the gauge invariance of the zero 
mode has a special nature. As we will show, these further terms are responsible for current 
renormalization. To simplify the discussion, we will restrict ourselves to the effective theory 
operators for zero frequency photons. 

Using gauge symmetry, we can write down the general form of the ultra low-energy 
effective theory for a single (p field coupled to zero frequency photons. We choose to construct 
this theory using Wilson lines. By cycling once over the z-th compact dimension, we can 
form gauge invariant Wilson lines W; of the form, 

Wi = exp j> dx % AjJ . (22) 

Notice that there is no sum over repeated indices in this definition. Due to the periodicity 
of the gauge field in the i-th direction, the loop integral 

/ dxiAi(x,t) = LAi(t) +L M^ty mn ' x/L , (23) 

produces just the rij = modes of the gauge field. Indeed the gauge transformation of the 
zero and non-zero modes, Eq. f[T7|) . demonstrates that the Wilson line Wj is gauge invariant. 
For our purpose, we wish to isolate completely the gauge field zero-mode and accordingly 
form modified Wilson lines Wi given by 

Wi = VoWiVl (24) 

where Vq is an operator that projects onto the zero-mode of the gauge field. A practical 
way to implement the action of Vo is to change the loop integration 

dxi — > <J> \\dx,j = — j> dx (25) 



so that 



Wi = exp (i j> dx Ai/L 2 ^ . (26) 

Furthermore it is useful to define Hermitian combinations of modified Wilson lines that 
transform simply under parity and charge conjugation, 

M +) = \(Wi + W}) (27) 



Notice that because W^ — y 1 — [W} ^] 2 , any operator involving W^ can be traded in 

for a tower of operators involving w/ . Hence we can build our theory solely in terms of 
W^ operators. 
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In addition to gauge, C, P, and T invariance, the theory on a torus has S4 cubic invari- 
ance. Writing down operators consistent with these symmetries, we arrive at the following 
ultra low-energy effective Lagrangian for a single ip field 



m 



|2 



C = \D^\ 2 -m(L) 2 |yf + Q 1 (L)W^ ) J + Q 2 {L) (W (_) • W H ) 

+ Q 3 (L) (WH ■ Dip*) (W(-) • Dcp) + Q 4 (L) J^W^W^ 'A^A^ 

+ Q 5 (L) (W<~> ■ W^) (W H • J) + Q 6 (L) ^ V^ ( ~ ) V^ H W i ( " ) J< + . . . . (29) 

i 

Above we have employed the current operator J, given by J = i [(Dip*)ip — ip*(Dip)]. A 
number of things about Eq. (1291) must be clarified. The . . . denotes that we have not 
finished writing the general Lagrangian allowed by symmetries. The most general Lagrangian 
contains a tower of terms with n insertions of operators. Writing down all such terms 
consistent with S4 for a given n is arbitrarily complicated. Fortunately the series expansion of 
w\ ■* in terms of the gauge field starts out at a single zero-frequency photon. Thus operators 
with n insertions of w\ ^ contribute to processes with at least n zero-frequency photons. In 
Eq. (j29p . we have written down all operators with at most three insertions of W*>~\ Thus 
the Lagrangian generates all possible couplings to at most three zero-frequency photons. We 
have also restricted the dynamics to ultra-low energies, so have only kept terms with at most 
two derivatives, D, acting on ip. Finally while a term of the form, (W^ ■ W^) iD^ip] 2 , is 
allowed by symmetries, it has been removed by a field redefinition. 

The coefficients Qj{L) in Eq. (|29|) must be determined from matching, and thus in general 
require the calculation of loop graphs with an arbitrary number of photons in the microscopic 
theory. It is possible that certain additional symmetries of the underlying theory constrain 
some coefficients to vanish. Because this is the zero-frequency sector of an effective theory 
for a stable particle, no multi-particle production thresholds can be attained in loop graphs 
that determine the matching coefficients. Thus in a sympto tically large volumes, the new 



coupling constants Qj(L) will be exponentially small [271. 1281. l29j|. Consequently SO (4) will 
be restored in large volumes. 



C. Zero frequency effective theories 

Using the general analysis from above, it is straightforward to construct single particle 
effective theories that reproduce the zero-frequency results derived in Sec. [TTJ There is one 
difference, however. The underlying theory, QCD, has quark fields with fractional charges. 
Maintaining periodicity of both quark fields under zero-mode gauge transformations requires 
a slightly modified quantization condition, namely 

67T 

ex = —n. (30) 
eL 

This modification reflects that both quark charges are quantized in units of e/3. The Wilson 
lines Wj are now given by 

W l = exp I j & dxi Ai J , (31) 
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and similarly for the W\. Thus for charged and neutral pions, 6 we require the effective 
Lagrangian 

C = -tr(D^Dv$) - -m 7r (L) 2 tr($ 2 ) + iQ(L) ■ tr [Q(D$)$ - Q$(D$)] 

2 2 



+ Q{L) 2 W { - ) ■ W { -hr 



(g$) 2 + V$ 2 

5 



Q{Lf w H ■ w { - ] tr [(g$) 2 - g 2 $ 2 ] , 

(32) 



where m n (L) includes the infrared running of the pion mass (not calculated here), and the 
new coupling constants Q(L), Q(L), and Q(L) are given by 

Q(L) = -7^1 [2/i/2(m 2 , L) + mlh /2 {ml L)] (33) 

Q(L) 2 = ^KhMmlL) (34) 

Q(L) = 0. (35) 

For the charged pions, we have also calculated all two-photon graphs to one-loop order using 
Eq. (j2j) (this includes the one-pion irreducible contributions shown in Fig. [31 and additionally 
the set of one-pion reducible diagrams which are not depicted) and do not find the need for 
an extra two-photon coupling to charged pions in Eq. (1321) . For this reason, the coupling 
constant Q(L) vanishes. The term with coefficient Q(L) only couples zero-frequency photons 
to neutral pions. 

The single particle effective theory described by Eq. (|32j) correctly reproduces the in- 
frared running of one- and two-photon processes for both charged and neutral pions at zero 
frequency. This theory is gauge invariant in finite volume because of the allowance for new 
operators which are Wilson lines that cycle the compact dimensions. These operators, more- 
over, lead to violation of 5*0(4) invariance. Consider a charged pion at rest, P = 0. The 
current in Eq. (jHJ) is 

(7^(0)1^1^(0)) = 2m 7T {±e)g fl0 , (36) 

where 2m n is a relativistic normalization factor. Boosting to a frame where P^0, generates 
a current 

r±(r>\\ ]T|J:/P^ — O. /P2 i m 2 



**(P)\ J|tt ± (P)) = 2VP2 + ml (±e)[l - Q(L)} V, (37) 



where V = P '/ 'y P 2 + m 2 is the relativistic velocity. Because of 5*0(4) breaking, the current 
in this frame is not simply the charge times the velocity, J ^ (±e)V. Instead, the current 
is screened by finite volume effects, J = (±e)[l — Q(L)]V. 



IV. FIELD THEORY IDENTITIES 

Above we have derived finite volume modifications to the current of the charged pions 
and the zero-frequency scattering tensor for the neutral pions. While we have accounted for 



6 Because electromagnetism explicitly breaks isospin symmetry, we should formulate the low-energy theories 
for charged and neutral pions separately. Although we utilize traces over SU{2) pion fields Q(x,t), there 
are no interactions between pions in Eq. (|32p. 
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these findings using gauge invariant single particle effective theories, here we show that our 
results are completely consistent with field theoretic identities valid in finite volume. 



A. Electromagnetic Vertex 

The zero-frequency part of the electromagnetic vertex is constrained by gauge invariance 
via the Ward identity [30j. Let r M (P, P) denote the zero- frequency electromagnetic vertex 
function of the charged pions. The Ward identity requires 

- t T»(P,P) = Qe^-G(P)-\ (38) 

where G(P) is the pion propagator. In finite volume, we found that the wave function 
correction did not exactly cancel the forward part of the vertex function. This lead to the 
new coupling Q(L) in Eq. (1331 . Thus at finite volume, the differential form of the Ward 
identity shown in Eq. (1381 is violated. Quite simply, however, the steps used to derive 
Eq. (|38|) are not valid in a fixed finite volume. 



On the other hand, starting from the Ward-Takahashi identity [31|, [32J we have 



- ik^{P + k,P) = Qe [G(P + A;)" 1 - G(P)- 1 ] . (39) 

This identity is valid in finite volume. We can demonstrate this explicitly using the charged 
pion vertex function, r M (P + A;, P) = (tt ± (P + k)\ J M |7r ± (P)). To one-loop order, we evaluate 
the diagrams in Fig. [1] and contract with the momentum transfer, k^. We find 



-a>(7r ± (P + A;)|J^|7r ± (P)) 
= =pe|(2P + k) ■ k - 

where we have abbreviated 




(2P + k)-k (2P + k) ■ q (2q + k) ■ k 



q 2 — m\ [q 2 — ml] [(q + k) 2 — m 2 ] 



,(40) 



H q=2ivn/L 

and implicitly regulate ultraviolet divergences using dimensional regularization. We then 
write 

(2q + k) ■ k = (q + k) 2 - ml - [q 2 - m% 

in order to reduce factors in the numerator of the last term. To arrive at the Ward-Takahashi 
identity from Eq. (j40p . we must show that the terms in the J^vanish. This follows imme- 
diately by using discrete translational invariance to re-index the sum. If k ^ 2irm/L, then 
the summation over spatial momentum modes cannot be re-indexed in this manner. Con- 
sequently the validity of the Ward-Takahashi identity, Eq. (1351 . in finite volume hinges on 
quantized photon momentum. 

Having established that the Ward-Takahashi identity holds in finite volume, there must 
be a flaw in the subsequent derivation of the Ward identity. To arrive at the differential 
form of the identity, Eq. (|38|) . from Eq. (|39|) a limiting process k^ — > is required. At fixed 
volume, the spatial momentum quantization condition invalidates this procedure. Contrary 
to Eq. (!38l) . there is no condition imposed on r M (P, P) in a compact space. In finite volume 
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with infinite time extent, only the spatial part of the differential form of the Ward identity 
does not hold. One can take the limiting procedure with respect to the zeroth component of 
momentum transfer, k° — > 0. Consequently the time component of Eq. (!38|) remains valid. 
Our results are indeed consistent with this fact, cf. Eq. ([7]). 7 



B. Compton Tensor 

The classical Thomson cross section arises in the zero-frequency limit of electromagnetic 
waves scattering off charged particles. According to low-energy theorems 35, 36[, any sen- 



sible gauge invariant field theory of charged particles will reproduce the Thomson cross 
section. In terms of the off-shell Compton scattering amplitude for a scalar particle, the 
zero-frequency part is required to be of the form 

T, v {u = J = 0) = 2{Qef (g, v - ^f^j , (42) 

where P^ is the particle's four momentum. Upon squaring and multiplying with phase space 
factors, the first term produces the Thomson cross section, while the second term is the Born 
contribution (which survives when we take the zero frequency limit before going on-shell). 
For the neutral pion, Eq. (j4"2"l) mandates that the Compton tensor vanishes, contrary to our 
results in finite volume, Eq. (fTTj) . 

The Thomson limit of the Compton tensor can be derived rigorously in field theory from 
generalized Ward identities, specifically for a scalar particle we have 



iT^u = a/ = 0) = {QefG(P)- 1 



d 



.2 



dP^dP 



-MP) 



G(PY\ (43) 



which reproduces both the Thomson and Born terms. This generalized Ward identity for 
the two-photon amplitude is not valid in finite volume; because, as with its counterpart in 
Eq. (!38|) . its derivation relies on a limiting procedure. 

Returning to the step in the derivation of Eq. (l4"3l) before the limiting procedure, we have 
a version of the Ward-Takahashi identity that is valid in finite volume. Let the initial particle 
(photon) momentum be denoted by P (k), and the final particle (photon) momentum by P' 
(k'). Then we have 

k w P G(P')iT^(P', k'; k, P) G(P) = (Qef \G{P + k) - G(P') - G(P) + G(P - fc')l • (44) 

Using the analytic expression for the one-loop diagrams in Fig. [3] for the neutral pion, 
one can verify explicitly that Eq. (jUJ) holds in finite volume provided that the photon 
momenta, k and k', are quantized. The validity of Ward-Takahashi identities requires 
discrete translational invariance. 

Now by taking the limit k' — ► k, followed by k — > in Eq. (|44|) . we accordingly recover 
the differential form of the identity in Eq. (|43p . Quite simply then, the zero frequency part 



7 xPT studies of the volume effects for form factors of pseudoscalar mesons 33|, [3J] have utilized only 
the time-component of the current, and considered the extent of the time direction as infinite. In this 
framework, no modification to meson charges was found, consistent with Eq. 
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of the Compton tensor is not constrained in finite volume as the limit k — > cannot be 
taken. Gauge symmetry constrains only the frequency dependent combination appearing in 
Eq. (I44"j) . Because we have kept the time direction infinite, a limiting procedure does exist 
for the time-time component of the scattering tensor. Consequently Eq. (|43p must apply to 
T qo (uj = uj' = 0), as is indeed the case for our one-loop results for the neutral pion, Eq. (fiT|) . 8 

V. CONCLUSION 

Above we have considered infrared effects on currents in finite volume field theories. 
Using the chiral Lagrangian as an example, we showed that matter fields can additively 
renormalize electric current in finite volume. Such effects do not violate gauge invariance; 
on the contrary, new couplings are allowed because of periodicity constraints on zero-mode 
gauge field transformations. Consequently gauge invariant single particle effective theories 
can be formulated that reproduce the infrared behavior of the interacting theory. These 
theories are written in terms of Wilson lines that cycle over the compact dimensions. As 
SO (4) is explicitly broken in these theories, boosting a charged particle from its rest frame 
to a frame moving with velocity V does not result in a current J = QeV . There are no 
contradictions with Ward-Takahashi identities, or low-energy theorems. Differential forms 
of Ward identities are inapplicable in finite volume. 

Conserved currents are not protected from infrared renormalization in finite spaces with 
discrete translational invariance. As non-perturbative field theories, such as QCD, are nu- 
merically simulated in a finite Euclidean space, it is important to understand the infrared 
running of current couplings. As a practical application of our work, the single particle 
effective theory derived here can be extended to describe volume effects for properties of 
hadrons determined from lattice QCD. 
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By taking the limit k° — > in the singly contracted identity 

WG{P')iT^G{P) = Qe [G{P')iT y {P',P + k)G(P + k) - G{P' - k)iT v {P' - k, P)G(P)} , (45) 

we additionally see that Tq v (uj = u)' = 0) = for the neutral pion in finite volume with infinite time 
extent. Similarly the other singly contracted identity yields T^iui = w' = 0) = upon taking the limit 
k'° to zero. Both of these conditions are satisfied by our one- loop results, Eq. (jTTJ) . 
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FINITE VOLUME FUNCTIONS 



For processes without momentum insertion, all finite volume matching terms, AX(L), in 
Eq. ([6]) can be cast in terms of the basic building block 



Ip{m\L) = — J]— — - 



L 3 ^ [q 2 + m 2 f J (2vr) 3 [q 2 + 



dq 



1 



m 



21/3 



[m 



2\3/2-/3 poo 



(4vr)3/2r(/3) J 



(46) 



(47) 



where "d^{z, q) is a Jacobi theta function. To see that all other required finite volume func- 
tions can be written in terms of 7g(m 2 , L), we first define 



1 q 11 ■ ■ ■ q tn 



L 3 ' [qr 2 + m : 



21/3 



q %1 ■ ■ ■ q ln 
(2vr) 3 [q 2 + m 



21/3 • 



(48) 



As a consequence of cubic invariance in the sums, we have I % } "' l ™(m 2 , L) = 0, for n odd. For 
even values, we find 



/J- i2 "(m 2 ,L) 



2n 



L_ S {ni* . . . 5 nn-in n } ^ Q (-m 2 )^ +i _ n (m 2 , L). 



(49) 



The bracketed indices denote complete symmetrization in the usual way, e.g., {AB} 
UAB + BA). 
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